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The first measurement of direct photons in Au+Au collisions at ^/s NN = 200 GeV is presented. 
The direct photon signal is extracted as a function of the Au+Au collision centrality and compared 
to NLO pQCD calculations. The direct photon yield is shown to scale with the number of nucleon- 
nucleon collisions for all centralities. 

PACS numbers: 25.75.Dw 



One of the most exciting observations from experi- 
ments at the Relativistic Heavy Ion Collider (RHIC) is 
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the strong suppression of the yield of hadrons at large 
transverse momenta (px > 2 GeV/c) in central Au+Au 
collisions, as compared to measured yields in p + p colli- 
sions scaled by the number of binary nucleon-nucleon col- 
lisions [ESS E|- Such quenching was predicted to result 
from the energy loss of hard-scattered partons propagat- 
ing through the high density matter created in heavy ion 
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collisions [5j. It was later proposed that the observed 
hadron suppression could be an initial-state effect due 
to saturation of the initial parton distributions in large 
nuclei 0. The high-py hadron suppression was not ob- 
served in d+ Au collisions |7J, |8j . This indicates that the 
suppression in Au+Au collisions is due to the extended 
dense matter in the final state, that is absent in d+Au 
collisions. 

Measurement of direct photon production allows more 
definitive discrimination between initial- and final-state 
suppression due to the fact that photons, once produced, 
are essentially unaffected by the surrounding matter. 
Hence photons produced directly in initial parton scat- 
terings will not be quenched unless the initial parton dis- 
tributions are suppressed in the nucleus. In fact, there 
may be additional direct photon yield in AA collisions Q 
due to various processes such as momentum broadening 
of the incoming partons, additional fragmentation con- 
tributions jlOl Hll | , or additional scatterings in the ther- 
malizing dense matter of the final state. 

This letter reports on direct photon production in 
Au+Au collisions at ^Js NN = 200 GeV with data taken 
by the PHENIX experiment U3 during the second RHIC 
run (2001). This analysis used the Beam-Beam Counters 
(BBC, 3.0 < \r]\ < 3.9) and the Zero Degree Calorimeter 
(ZDC) for trigger and event characterization, the Electro- 
magnetic Calorimeter (EMCal) in the two central arms 
(|ry| < 0.35) to measure the inclusive 7, 7r°, and ?/ yields, 
and the tracking system of the central arms to estimate 
the charged particle contamination. The EMCal con- 
sists of two subsystems: six sectors of lead-scintillator 
sandwich calorimeter (PbSc) and two sectors of lead-glass 
Cherenkov calorimeter (PbGl). Located at a radial dis- 
tance of about 5 m each sector covers an azimuthal inter- 
val of A0«22.5°. The fine segmentation of the EMCal 
(AcfixArj ~ 0.01x0.01) ensures that the two photons from 
a decayed ir° are well-resolved up to transverse momenta 
of 15-20 GeV/c. 

The event centrality was selected by cuts on the cor- 
related distribution of charged particles detected in the 
BBCs versus energy measured in the ZDC detectors. A 
Glauber model Monte Carlo combined with a simulation 
of the BBC and ZDC responses gave an estimate of the 
associated number of binary collisions (N co u) and partic- 
ipating nucleons (N part ) for each centrality bin (values 
tabulated in Ref. Q)- 

For this analysis a mimimum bias trigger sample of 
30 x 10 6 events, also used for the previously published n° 
analysis was combined with a Level-2 trigger event 
sample equivalent to an additional 55 x 10 6 minimum 
bias events. The Level-2 trigger sample was obtained 
by use of an EMCal software trigger on highly energetic 
showers equivalent to the Level-1 hardware trigger used 
in Ref. [Lj. The threshold energy of the trigger was set 
at 3.5 GeV with a resulting trigger efficiency plateau at 
100% for single photons above pr ~ 5 GeV/c (6.5 GeV/c) 
for the PbSc (PbGl). The normalization of the Level-2 
data sample relative to the minimum bias data sample 



is accurate to 2%. In the following, the minimum bias 
result refers to the combined Level-2 and minimum bias 
trigger samples without selection on centrality. 

The direct photon yield is extracted on a statistical 
basis, without isolation cuts, by a comparison of the in- 
clusive photon spectra to the expected background from 
hadronic decays 0,^3 (mainly 7r° — * 27). Photon- like 
clusters are identified in the EMCal by applying appro- 
priate Particle Identification (PID) cuts based on time- 
of-flight and the shower profile. The consistency of the 
final results obtained independently with the PbSc and 
PbGl, and with different PID cuts, including no PID cut, 
is used to check the systematic error estimates. The 7r° 
and 77 yields are determined as described in 0, IT^ | by an 
invariant mass analysis of photon pairs, with the combi- 
natorial background established by combining uncorre- 
cted photon pairs from different events. 

The raw inclusive photon-candidate spectra must be 
corrected for charged and neutral hadron contaminations 
not removed by the PID cuts, as well as for photon con- 
versions. Charged contaminants are identified by associ- 
ating photon candidates in the EMCal with charged hits 
in the pad chamber (PC3) positioned directly in front of 
the EMCal. The charged contaminant spectra are sub- 
tracted from the photon-candidate spectra. The charged 
hadron contamination depends strongly on the PID cut 
and increases significantly for px < 3 GeV/c with a con- 
tribution of 4% above 3 GeV/c for the tightest PID 
cut. The contamination of neutral hadrons (mainly anti- 
neutrons) is determined with a full GEANT simulation of 
the detector response to neutrons and anti-neutrons with 
input spectra based on the proton and anti-proton yields 
measured by PHENIX |4|. The neutral hadron contam- 
ination is found to be negligible above pr — 5 GeV/c 
(< 1%). The neutral photon-candidate spectra are cor- 
rected for conversions removed by the charged contami- 
nant subtraction with a pr-independent factor (5.9-7.3% 
for different sectors based on simulation). 

The raw spectra are normalized to one unit of rapidity 
and full azimuth (the purely geometrical acceptance cor- 
rection is ~ 1/0.35) . The spectra are further corrected 
for (i) the detector response (energy resolution, dead ar- 
eas), (ii) the reconstruction efficiency (PID cuts), and 
(iii) occupancy effects (cluster overlaps). These correc- 
tions are quantified by embedding simulated single 7's, 
7r 0, s, or 77's from a full PHENIX GEANT simulation into 
real events, and analyzing the merged events with the 
same analysis cuts used to obtain the real yields. The 
overall ir° yield correction was ~2.5 with a centrality de- 
pendence of <25%. The losses were dominated by fidu- 
cial and asymmetry cuts. The nominal energy resolution 
was adjusted in the simulation by smearing the energies 
with a constant term of ~ 5% for PbSc and ~ 7% for 
PbGl to reproduce the measured width of the ir° peak 
observed at each px- The shape, position, and width of 
the 7T° peak measured for all centralities were confirmed 
to be well reproduced by the embedded data. 

The energy calibration of the EMCal was corroborated 



TABLE I: Summary of the dominant sources of systematic er- 
rors on the 7T° and inclusive 7 yields extracted independently 
with the PbGl and PbSc electromagnetic calorimeters. The 
error estimates are quoted at two pr values in central events 
for the PbGl and PbSc. For the combined 7r° and inclusive 7 
spectra and 7/71° ratios, the approximate statistical and sys- 
tematical errors are quoted for the most peripheral and most 
central reactions. 



7T° error source 


PbGl (Central) 
3GcV/c 8.5GeV/c 


PbSc (Central) 
3GeV/c 8.5GeV/c 


Yield extraction 
Yield correction 
Energy scale 


8.7% 
12.1% 
13.8% 


7.0% 
12.0% 
14.1% 


9.8% 
10.3% 
10.5% 


7.2% 
12.5% 
11.4% 


Total systematic 


20.3% 


19.8% 


17.7% 


18.4% 


Statistical 


10.6% 


32.5% 


2.1% 


10.5% 


Inclusive 7 error 


Non-7 correction 
Yield correction 
Energy scale 


2.4% 
10.2% 
15.7% 


2.4% 
12.0% 
13.7% 


3.2% 
9.1% 
12.4% 


3.2% 
12.5% 
10.8% 


Total systematic 


18.9% 


18.4% 


15.7% 


16.8% 


Statistical 


1.2% 


14.1% 


0.6% 


4.1% 


7/jr syst. 
7/7T stat. 


13.6% 
10.7% 


12.6% 
35.4% 


14.0% 
2.2% 


13.4% 
11.3% 


Total errors PbGl and PbSc 


combined 




Error 


Peripheral 
3GeV/c 8.5 GeV/c 


Central 
3GeV/c 8.5GeV/c 


7T° syst. 
7T° stat. 


13.2% 
3.0% 


17.0% 
35.3% 


13.9% 
1.8% 


16.1% 
9.6% 


7 syst. 
7 stat. 


11.4% 
3.0% 


15.6% 
28.8% 


11.5% 
0.6% 


15.9% 
3.8% 


7/tt° syst. 
7/-7T stat. 


9.9% 
4.2% 


13.1% 
45.6 % 


9.7% 
1.9% 


11.2% 
10.3% 


7/7T bkg calc. 




4% 




4% 



by the position of the 7r° invariant mass peak, by the 
energy deposit from minimum ionizing charged particles 
traversing the EMCal (PbSc) , and by the correlation be- 
tween the measured momentum of electron and positron 
tracks identified by the ring-imaging Cherenkov detector 
and the associated energy deposit in the EMCal. From 
these studies it is determined that the accuracy of the 
energy measurement was better than 1.5%. 

The main sources of systematic errors in the PbSc and 
PbGl measurements are the uncertainties in: (i) the yield 
extraction, (ii) the yield correction, and (iii) the energy 
scale. The relative contributions of these effects to the 
total error differ for the PbSc and PbGl (Tabled. The 
weighted average of the two independent measurements 
reduces the total error. The final systematic errors on 
the spectra are at the level of ~ 15 — 20% (Table 0. A 
correction for the true mean value of the px bin is applied 
to the steeply falling spectra. 

The completely corrected and combined PbSc and 
PbGl inclusive photon yields are compared to the ex- 
pected yields of background photons from hadronic de- 
cays in Fig. n f° r minimum bias Au+Au collisions and 
for five centrality bins. The decay photon calculations 
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FIG. 1: Double ratio of measured (7/71"°) Measured invariant 
yield ratio to the background decay (7/7r°)Background ratio as 
a function of pr for minimum bias and for five centralities of 
Au+Au collisions at ^/s~^ = 200 GeV (0-10% is the most 
central). Statistical and total errors are indicated separately 
on each data point by the vertical bar and shaded region, re- 
spectively. The solid curves are the ratio of pQCD predictions 
described in the text to the background photon invariant yield 
based on the measured n° yield for each centrality class. The 
shaded region around the curves indicate the variation of the 
pQCD calculation for scale changes from p-r/2 to 2px, plus 
the (Ncoll) uncertainty. 



are based on the measured 7r° and rj spectra assum- 
ing mr-scaling for all other radiative decays (r]',Kg,ui). 
The comparison is made as the ratio of measured (in- 
clusive) 7/71"° and calculated background 7/77° since this 
has the advantage that many uncertainties, such as the 
energy scale, cancel to varying extent in the ratio. Since 
the 7T° spectra of the background calculations are taken 
to be the same as the measured spectra we have 



(7/77°) 

D V [ ' / Measured 

7 " (7 A ) 



TMcasurcd 



Background ^Background 



(i) 



and any significant deviation of the double ratio above 
unity indicates a direct photon excess. In Fig.^an excess 
is observed at high px with a magnitude that increases 
with increasing centrality of the collision. 

The measured results are compared to NLO pQCD 
predictions [Tt^ , scaled by the number of binary nucleon 
collisions for each centrality selection. The same calcula- 
tions are in a gree ment with the PHENIX direct photon 
measurement [l5j for p+p collisions at the same y/s, and 
similar NLO pQCD calculations provide a good descrip- 
tion of the measured 77° production in p+p collisions [l3| . 
The calculations were performed |l5l Il7j with normaliza- 
tion and factorization scales equal to px, the CTEQ6 0] 
set of parton distribution functions, and the GRV set of 
fragmentation functions 19]. The direct photon spectra 
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FIG. 3: Ratio of Au+Au yield to p+p yield normalized by the 
number of binary nucleon collisions as a function of centrality 
given by N par t for direct 7 (closed circles) and n° (open cir- 
cles) yields integrated above 6 GeV /c. The p+p direct photon 
yield is taken as the NLO pQCD prediction described in the 
text. The error bars indicate the total error excluding the 
error on (N co u) shown by the dashed lines and the scale un- 
certainty of the NLO calculation shown by the shaded region 
at the right. 



FIG. 2: Direct 7 invariant yields as a function of trans- 
verse momentum for 9 centrality selections and minimum bias 
Au+Au collisions at ^s NN = 200 GeV. The vertical error 
bar on each point indicates the total error. Arrows indicate 
measurements consistent with zero yield with the tail of the 
arrow indicating the 90% confidence level upper limit. The 
solid curves are pQCD predictions described in the text. 



extracted as 7 D iroct = (1 - R y x ) • TMeasured are shown in 
Fig. [2 for all nine centrality selections as well as mini- 
mum bias, and compared to the same NLO calculations. 
The binary collision scaled predictions are seen to provide 
a good description of the measured direct photon spec- 
tra (Fig. 0). The increasing ratio with centrality seen 
in Fig. is therefore attributed to the decreasing decay 
background due to it suppression Q. 

Medium effects in AA collisions are often presented 
using the nuclear modification factor given as the ratio 
of the measured AA invariant yields to the TV TV-collision- 
scaled p + p invariant yields: 



Raa(pt) = 



(l/N e /i)d 2 N AA /d PT dy 



{N coU }/a p nel 



vi> 



d 2 <r pp /dp T dy' 



(2) 



where the (N co u) / a p ™ el is the average nuclear thickness 
function, (Taa), in the centrality bin under consideration 
(Ref [3]). Raa(pt) measures the deviation of A A data 
from an incoherent superposition of NN collisions. 

The centrality dependence of the high px 7 production 
represented as a function of the number of participating 
nucleons, N part , is shown by the closed circles in Fig.|3] 
The production in Au+Au collisions relative to p + p is 
characterized by the Raa(pt > 6 GeV/c) ratio of Eq. @ 
as the ratio of Au+Au over the (N co u) -scaled p+p yields 



each integrated above 6 GeV/c. The direct photon p + p 
yields are taken as the NLO pQCD predictions described 
above. As noted above, the high px direct 7 production 
is observed to scale as the (N co u) -scaled p+p 7 yieldpre- 
diction for all centralities. This is in sharp contrast 0] to 
the centrality dependence of the ir° Raa(pt > 6 GeV/c) 
shown by open circles in Fig. [31 where the measured ir° 
yield [13| is used as the p + p reference in Eq. @ . 

The observed close agreement between the measured 
yields and NLO calculations is in striking contrast to 
observations for central Pb+Pb collisions at -Js NN = 
17.3 GeV 01 where the measured photon yield exceeds 
the (iV co /;)-scaled p + p yield by about a factor of two. 
The present result constrains modifications of the initial 
parton distributions, or of the fragmentation contribu- 
tions [Tol HH (in these NLO calculations the contribu- 
tion is significant: ~ 50% at 3.5 GeV/c and ~ 35% at 
10 GeV/c), or additional photon yield from thermal ra- 
diation to levels comparable to the present measurement 
uncertainty. 

In summary, the transverse momentum spectra of di- 
rect photons have been measured at mid-rapidity up to 
Pt ~ 13 GeV/c for nine centrality bins of Au+Au colli- 
sions at ^ys NN = 200 GeV. The significance of the direct 
photon signal increases with collision centrality due to 
the increasingly suppressed ir° production and associated 
decrease in the photon background from hadron decays. 
The direct photon spectral shapes and invariant yields 
are consistent with NLO pQCD predictions for p + p re- 
actions scaled by the average number of inelastic NN 
collisions for each centrality class. The close agreement 
between measurement and the binary scaled pQCD pre- 
dictions of the direct photon yield suggests that nuclear 
modifications of the quark and gluon distribution func- 
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tions in the relevant region of momentum fraction x are 
minor. The result provides strong confirmation that the 
observed large suppression of high px hadron production 
in central Au+Au collisions is dominantly a final-state 
effect due to parton energy loss in the dense produced 
medium, rather than an initial-state effect. 
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